function T = ECI2SCI(year)
%[] Using the given year, calculate the transformation matrix from ECI to
%SCI.

obliquity_dms = [23 , 27 , 8.26-(0.4684*(year-1900))/216000];
%[degrees minute seconds] obliquity of ecliptic

obliquity_deg = Dms2Deg(obliquity_dms)*pi/180;
%[deg] Obliquity of ecliptic in degrees.

T = [1,0,0;...
     0,cos(obliquity_deg),sin(obliquity_deg);...
     0,-sin(obliquity_deg),cos(obliquity_deg)];
%[] Transformation matrix from ECI 2 SCI
end
